Counterpropagating optical beams and solitons 
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Physics of counterpropagating optical beams and spatial optical solitons is reviewed, including the 
formation of stationary states and spatiotemporal instabilities. First, several models describing the 
evolution and interactions between optical beams and spatial solitons are discussed, that propagate 
in opposite directions in nonlinear media. It is shown that coherent collisions between counterpropa- 
gating beams give rise to an interesting focusing mechanism resulting from the interference between 
the beams, and that interactions between such beams are insensitive to the relative phase between 
them. Second, recent experimental observations of the counterpropagation effects and instabilities 
in waveguides and bulk geometries, as well as in one- and two-dimensional photonic lattices, are 
discussed. A variety of different generalizations of this concept are summarized, including the coun- 
terpropagating beams of complex structures, such as multipole beams and optical vortices, as well 
as the beams in different media, such as photorefractive materials and liquid crystals. 

PACS numbers: 42.65.Jx, 42.65.Tg, 42.65.Sf, 42.70.Mp. 



I. INTRODUCTION 

One of the simplest processes in nonlinear (NL) optics 
leading to a variety of complex NL physics phenomena is 
the mutual interaction of two counterpropagating (CP) 
optical beams in a NL medium, capable of nonlinearly 
changing the refractive index of the medium. The under- 
lying geometry is conceptually very simple (see Fig. [1]): 
two beams enter a finite NL medium from the opposite 
sides and, when they overlap by their evanescent fields, 
the beams start interacting via the mutual NL change 
of the optical refractive index. A configuration of two 
waves interacting in a NL material is one of the most fre- 
quent ones in laser physics and wave mixing experiments. 
Numerous concepts in NL optics, such as phase conjuga- 
tion, Bragg reflection by volume gratings, wave-mixing 
in photorefractives, etc., are based on this simple geom- 
etry. Nevertheless, this simple geometry can give rise 
to an extremely complicated and sometimes counterin- 
tuitive dynamical behavior, including both mutual beam 
self-trapping and the formation of stationary states, as 
well as complex spatiotemporal (ST) instabilities It 
is for these reasons that CP beam configurations have 
achieved a paradigmatic role in NL physics of optical 
systems. 

Instabilities and chaos are typically expected to appear 
in NL optical systems that feature coupling as well as 
feedback as necessary ingredients. Therefore, CP waves 
were first studied in more complex systems than the con- 
figuration described above, foremost in NL optical res- 
onators. CP beams in Fabry- Perot resonators have been 
studied by Ikeda et al. 0, @l • They demonstrated that 
a ring resonator with Kerr nonlinearity undergoes a se- 
ries of bifurcations, as the incident power is increased, 
leading to chaos and "optical turbulence". Similar ST 
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FIG. 1: General schematic of two counterpropagating waves 
interacting in a finite-length nonlinear medium. 



instabilities were observed later in simpler configurations 
with only a single mirror, known as the single feedback 
systems, especially in the ones exhibiting saturable non- 
linearities, such as atomic vapors 0, 0, [y, 0] > liquid crys- 
tals [8[, and photorefractive (PR) media [9(. For example, 
instabilities in the polarizations of CP beams were ob- 
served in atomic sodium vapor 

ill 

, when for higher in- 
tensities the polarization first varies periodically and then 
the system dynamics becomes chaotic (see also Ref. p~olj ) . 
Silberberg and Bar Joseph [l[ were the first to demon- 
strate that even without an external feedback, instabili- 
ties and chaos can be observed in the simplest geometry 
shown in Fig. [T] In their configuration, the origin of 
instability was the combined action of gain due to four- 
wave mixing and distributed feedback due to scattering 
from the grating formed by the interference of the inci- 
dent laser beams. This discovery opened the possibility 
of considering CP beam configuration as a fundamen- 
tal configuration to investigate NL physics phenomena 
in optics. 

Such a possibility was backed up by many subse- 
quent observations, showing that when taking into ac- 
count transverse extent of the interacting beams, an even 
more complicated behavior can emerge, including oscil- 
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latory transverse instabilities [ll|, ll2| and pattern for- 
mation [TH, [3 that require neither cavity nor fi- 
nite response time. Sometimes, these instabilities could 
be associated with interesting transverse structures, such 
as polarization domain walls (lrl [ijj , supported by the 
mixing of two CP laser beams in a NL isotropic dielec- 
tric medium. Since transverse instabilities of CP opti- 
cal beams lead to pattern generation, a natural question 
arises in CP systems: which patterns may survive in this 
geometry. Among the various pattern classes especially 
singulary structures and localized states gained interest, 
owing to their strong NL nature most often associated 
with subcriticial bifurcations [Lsl ]. 

In addition to these localized states or dissipative (cav- 
ity) solitons, which are stabilized due to the common ac- 
tion of gain, loss, nonlinearity and diffraction, the prop- 
agation and self-action of a single beam in a NL medium 
is also known to generate spatial solitons Therefore, 
a natural extension was to study the propagation and in- 
teraction of CP spatial solitons and their generalizations, 
and to build a connection between these solitons and the 
feedback or cavity solitons, thereby attempting to answer 
fundamental questions about the nature and relations of 
propagating and dissipative solitons in media with NL 
optical refractive index. 

Collisions between solitons are perhaps the most fasci- 
nating feature of soliton phenomena, because the inter- 
acting self-trapped wave packets exhibit many particle- 
like features [20| . Solitons that propagate in the opposite 
direction enable a natural mechanism of soliton collis- 
sions, resulting from the strong interaction of the two 
beams. The CP solitons interfere and give rise to an ef- 
fective grating. For copropagating solitons, the grating is 
periodic in the transverse direction, with a period much 
greater than the optical wavelength; thus the interacting 
solitons go through very few grating periods. For CP col- 
lisions in contrast the grating is in the propagation direc- 
tion, hence the interacting solitons go through many pe- 
riods. Consequently, the interaction in the CP scheme is 
strongly affected by the mutual Bragg scattering. More- 
over, in the case of incoherent interacting beams, the CP 
scheme also allows for strong interaction, due to cross- 
coupling of the beams via the common refractive index 
structure that is not present in the copropagating case. 

Owing to these reasons, CP solitons assumed a 
paradigmatic role in the physics of NL optics: mutual 
self-trapping of two CP optical beams was shown to lead 
to the formation of a novel type of vector (or bimodal) 
solitons pU . [22T ]. for both coherent and incoherent in- 
teractions. A more detailed analysis [23[ revealed that 
these CP solitons may display a variety of instabilities, 
accompanied by nontrivial temporal and spatial dynam- 
ics, leading to many subsequent studies devoted to this 
fundamentally new subject. 

This paper aims to review the fundamental physics of 
CP optical beams and spatial optical solitons, and their 
paradigmatic role in NL optics, including the whole range 
of NL physics phenomena, from the formation of station- 



ary states up to ST instabilities. It summarizes a number 
of recent important results for the evolution and interac- 
tion of optical beams and spatial solitons that propagate 
in opposite directions, thereby emphasizing their general 
importance for NL physics. 

The paper is organized as follows. In Sec. 2 we 
present the derivation of a one-dimensional (ID) model 
for the beam propagation in a planar structure in PR 
NL crystals, and then apply it to the analysis of mutual 
self-trapping and modulational instabilities (Mis) of CP 
beams and spatial solitons. Section 3 is devoted to the 
analysis of 2D models, where we discuss nontrivial rota- 
tional beam dynamics and the transverse pattern forma- 
tion. In Sec. 4 we present the key experimental results 
for both one- and two-dimensional geometries. The more 
special case of solitons counterpropagating in optical lat- 
tices is discussed in Sec. 5, where we summarize both 
theoretical and experimental results. Section 6 is devoted 
to the discussion of various generalizations of the concept 
of beam counterpropagation, including the counterprop- 
agation of multipoles and vortex optical beams, as well 
as the beam interaction in liquid crystals. Finally, Sec. 
7 concludes the paper. 



II. ONE-DIMENSIONAL SYSTEMS 

A. Theoretical models and background 

Early theoretical descriptions of CP self-trapped 
beams, in one transverse dimension and steady-state, 
were given in (2^|, where bimodal CP solitons in Kerr 
media have been treated, and in [2lj . where collisions 
of solitons propagating in opposite directions, in both 
Kerr and local PR media, have been addressed. Follow- 
ing a more general exposition (23j, we present here the 
basic equations for the propagation and interactions of 
CP beams in saturable PR media. The temporal behav- 
ior of CP self-trapped beams is included in the equations 
by a time-relaxation procedure for the formation of space 
charge field and refractive index modulation in the crys- 
tal. 

We consider two CP light beams in a PR crystal, in the 
paraxial approximation, under conditions suitable for the 
formation of screening solitons. The optical field is given 
as the sum of CP waves F exp(ikz + iu>t) + B exp(—ikz + 
iuit), k being the wave vector in the medium, F and B 
are the slowly varying envelopes of the beams. The light 
intensity / is measured in units of the background light 
intensity, also necessary for the generation of solitons. 
After averaging in time on the scale of the response time 
To of the PR crystal, the total intensity is given by 

1 + 1 = (1 + I a ) {1 + s [m exp(2ikz) + c.c] /2} , (1) 

where I = \F\ 2 + \B\ 2 , m = 2FB*/{\ + I ) is the mod- 
ulation depth, and c.c. stands for complex conjugation. 
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Here the parameter s measures the degree of temporal 
coherence of the beams relative to the crystal relaxation 
time. For e = 0, i.e., when the relative phase of the 
beams varies much faster than To, the beams are effec- 
tively incoherent. In the opposite case e = 1, the in- 
tensity distribution contains an interference term that is 
periodically modulated in the direction of propagation z, 
chosen to be perpendicular to the c axis of the crystal, 
which is also the x axis of the coordinate system. Beams 
are polarized in the x direction, and the external elec- 
tric field E e , necessary for the formation of self-trapped 
beams, also points in the x direction. The electric field 
in the crystal couples to the electrooptic tensor, giving 
rise to a change in the index of refraction of the form 
An = —n^r e ffE/2, where no is the unperturbed index, 
r e ff is the effective component of the electro-optic ten- 
sor, and E is the x component of the total electric field. 
It consists of the external field and the space charge field 
E sc generated in the crystal, E = E e + E sc . 

The intensity modulates the space charge field, which 
is represented in the normalized form 

E sc /E e = E + ^ [Ei exp(2ikz) + c.c] , (2) 

where Eq is the homogeneous part of the x component of 
the space charge field, and E\(x, z) is the slowly varying 
part of the space charge field, proportional to e. It is Eq 
that screens the external field, and E\ is the result of the 
interference pattern along the z direction. 

In the isotropic approach, one assumes a local approx- 
imation to the space charge field, and looks for a solution 
with the saturable nonlinearity E = E e /(1 + 1). Substi- 
tuting Eqs. (H|) and @ in this expression, and neglecting 
higher harmonics and terms quadratic in m, the steady- 
state solutions Eq = —Io/(l+Io) and E\ = — em/(l + / ) 
are obtained. The temporal evolution of the space charge 
field is introduced by assuming relaxation-type dynamics 

Td t E +E = -— - ( 3a ) 

rd t E x + E x =-^- , (3b) 
t + io 

where the relaxation time of the crystal r is inversely 
proportional to the total intensity r = tq/(1 + /), i.e., 
the illuminated regions in the crystal react faster. The 
assumed dynamics is that the space charge field builds 
up towards the steady state, which depends on the light 
distribution, which in turn is slaved to the slow change 
of the space charge field. As it will be seen later, this 
assumption does not preclude a more complicated dy- 
namical behavior. 

Selecting synchronous terms in the NL paraxial wave 
equation, leads to the propagation equations in the form: 

id z F + d 2 x F = T [E F + ExB/2] , (4a) 
-id z B + d 2 x B = T[E B + E{F/2] , (4b) 
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FIG. 2: (a)-(c) Coherent interactions between the CP (a),(c) 
forward and (b) backward solitons. (d)-(f ) Incoherent inter- 
actions between the CP (d),(f ) forward and (e) backward soli- 
tons. For comparison, interactions between coherent in-phase 
(g) and 7r out of phase (h), and incoherent (i) copropagating 
solitons. The plots show intensities. The arrow indicates the 
propagation direction of each beam. The propagation dis- 
tance is 2.5 Lo- Adopted from [2lj . 

where the parameter V = (knoXo) 2 r e f fE e is the dimen- 
sionless coupling strength, and the scalling x — > x/xq, 
z — * z/Ln, (F,B) — * (F, B) exp(— iTz) is used. Here x$ 
is the typical beam waist and Lo — 2kx\ is the diffrac- 
tion length. Propagation equations are solved numer- 
ically, concurrently with the temporal equations. The 
numerical procedure consists in solving Eqs. ^ for the 
components of the space charge field, with the light fields 
obtained at every step as guided modes of the induced 
common waveguide. 



B. Counterpropagating solitons 

We consider first the interactions in a configuration 
where the two CP beams are launched parallel to each 
other, but with a transverse spacing between them [2l| . 
The parameters are chosen such that the formation of 
spatial solitons is preferred. The coherent interaction 
between these parallel CP beams is shown in Figs. Uta)- 
HIc). Figures[2fd)-mf) show an incoherent interaction be- 
tween the same beams. For comparison, the same beams 
in copropagating scheme are simulated in Figs. [2Ig)-[2{i)- 
Figure [2fg) [DJh)] shows a coherent interaction in which 
the relative phase between the launched beams is [tt]. 
Figure [^i) shows an incoherent interaction. Clearly, the 
outcome of the interaction between the beams in the CP 
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FIG. 3: Bidirectional waveguide, (a) Total intensity distri- 
bution; (b) right-propagating and (c) left-propagating beams. 
Parameters: e = 1, T = 5, initial peak intensities |fb| 2 = 
\B L \ 2 = 1. The size of data windows is 10 beam diameters 
transversely by 2 diffraction lengths longitudinally. Reprinted 
from HI]. 



scheme is very different from the one in the copropagating 
scheme, in both the coherent and incoherent cases. First, 
in the copropagating scheme, the mutual force between 
the solitons is proportional to the relative phase between 
them, hence the interaction can be attractive [Fig. HJg)] 
or repulsive [Fig. [Hh)]. In contrast, in the CP case the 
relative phase oscillates on a scale much shorter than the 
soliton period, thus the relative phase does not play any 
role. 

The second major difference between the counter- and 
copropagating cases has to do with the radiation. The co- 
herent interaction in the CP scheme radiates [Figs. 
and[5Jb)], which again proves that this system is nonin- 
tegrable. On the other hand, the incoherent interaction 
between the CP solitons does not radiate much [Figs. 
[5Jd)-[Uf)]. Finally, one can see that a portion of the for- 
ward beam couples into the region where the backward 
beam is propagating. In the incoherent interaction, the 



forward beam gradually tunnels into the backward soli- 
ton region, hence the forward intensity at the backward 
soliton region increases monotonically [Fig. |2ff)]. This 
behavior represents an example of directional coupling 
or resonant tunneling. For coherent interactions the dy- 
namics are more complex, as the intensity coupled from 
the forward beam to the region "under" the backward 
beam oscillates [see the sidebands in Fig. H[c)], and, in 
contradistinction to the incoherent case, light does not 
accumulate in the "sidebands." 

Head-on collision of the beams with initial soliton pro- 
files, after temporal relaxation to a steady state, results in 
the formation of a CP soliton [23[ . Shooting initial beams 
with arbitrary parameters generally leads to the z depen- 
dent or nonstationary character of the beam propagation. 
In some domain of the initial parameters, for example 
with the relative angle of beam scattering 9 close to it and 
small initial transverse offset, the time-relaxation proce- 
dure converges to the stationary in time structures, which 
are identified as the steady-state self-trapped waveguides, 
or as bent CP solitons. The formation of a single bidi- 
rectional waveguide is shown in Fig. [31 Two coherent 
Gaussian beams are launched at different lateral posi- 
tions perpendicular to the crystal edges, 9 = tt. Both 
beams diffract initially, until the space charge field is de- 
veloped in time, to form the waveguide induced by the 
total light intensity, Fig. [3Ja). This induced waveguide 
traps both beams, Figs. [3jb) and[3^c). When the initial 
transverse separation is four or more beam diameters, the 
beams hardly feel the presence of each other, and focus 
into individual solitons. For the separation of two beam 
diameters, the interaction is strong enough for the beams 
to form a joint waveguiding structure, as is shown in Fig. 

El 



C. Splitup transition 

Consideration of a wider region of control parameters 
leads to a more complex picture. To capture the tran- 
sition from a CP soliton to a waveguide more clearly, 
the head-on collision of two identical Gaussian beams 
was considered in [H, HH ■ ^ n the absence of the other, 
each beam focuses into a soliton. The situation when 
they are both present, and when the coupling constant 
r and the crystal length L are both varied, is displayed 
in Fig. [U It is seen that in the plane (L,T) of control 
parameters there exists a critical curve below which the 
stable CP solitons exist (the first curve in Fig. At 
that critical curve a new type of solution appears, after a 
symmetry breaking transition, in which the two compo- 
nents no longer overlap, but split and cross each other. 
We term this phenomenon the splitup transition [23|. A 
few examples are depicted in the insets in Fig. 2) As 
the beams split, a portion of each beam remains guided 
by the other, forming bidirectional waveguides. Both the 
solitons and the waveguides are steady-state solutions. 
As one moves away from the first critical curve, into the 
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FIG. 4: Critical curves in the parameter plane for the ex- 
istence of stable CP solitons, bidirectional waveguides and 
unstable solutions. Below the first curve CP solitons exist, 
between the curves bidirectional waveguides appear. At and 
above the second curve unstable solutions emerge. Insets de- 
pict typical beam intensity distributions in the (a;, z) plane at 
the points indicated. The points are numerically determined, 
the curves are inverse power polynomial fits. Reprinted from 



id z F + -d 2 x F = Sn 2 F + -Sn 2 m B , (5a) 

-id z B+ 1 -dlB = §nlB+ 1 -{8nl i yF . (5b) 

To see the propagation behavior that is a mixture of the 
self-focusing and pattern formation, the counterpropa- 
gation of two wider beams in an anisotropic nonlocal 
medium is simulated in Fig. [5] Figures H^a) and [5Jb) 
show how the profiles of the beams change as they prop- 
agate. Figure [3(c) shows the profile of the forward beam 
as it leaves the crystal. It has split into three beams, rem- 
iniscent of the breaking of a uniform beam into stripes 
in the experiments on pattern formation in CP beams. 
The solid line in Fig. [DJd) shows the backward beam 
as it leaves the crystal. For comparison, the dashed line 
shows what the beam would look like if the nonlinearity 
were absent. One can see that on the one hand the back- 
ward beam gets amplified while propagating through the 
crystal; on the other hand the self-focusing effect of the 
nonlinearity is also clearly visible. The effect of the self- 
bending is weak, due to the short propagation distance. 
However, its effects are clearly visible in the asymmetry 
of the beam profile in Fig. Etc). 



region of high couplings and long crystals, a new critical 
curve is approached, where the steady-state waveguides 
loose stability. The second critical curve is also drawn in 
Fig. [4l and the insets to the curve show typical unsta- 
ble beam profiles. The shape of these curves suggests an 
inverse power law dependence, and the theory confirms 
such a dependence. At and beyond the second critical 
curve, dynamical solutions emerge. 



D. Anisotropic nonlocal theory 

Anisotropic nonlocal theory of the space charge field, 
induced by the coherent CP beams in biased PR crys- 
tals, is more involved than the isotropic theory (26|. It 
yields significantly different results from the isotropic lo- 
cal model, especially when the crystal c axis is tilted with 
respect to the direction of the propagation of the beams. 
A more complete description of CP beams requires inclu- 
sion of both the drift and the diffusion term. 

In the anisotropic approximation, the NL refractive 
index change 8n 2 can be decomposed into the form: 
8n 2 = Suq + 5ji 2 l [exp(2ikz) + cxp(— 2ikz)]/2. The mod- 
ulated (Sn^) and unmodulated (8n 2 ) parts are the func- 
tions of the space charge field and the nonzero compo- 
nents of the electro-optic tensor. The propagation equa- 
tions of the beam envelopes in the paraxial approxima- 
tion are now given by: 



E. Modulational instability 

As mentioned, the configuration of two waves inter- 
acting in a NL material is one of the most used in laser 
physics and NL wave mixing experiments. Instabilities, 
self-oscillations, and chaos, which are the fundamental 
processes in NL optics, can be observed very often in 
such systems. In ]1[ it was demonstrated, for the first 
time, that self-oscillations and chaos can be obtained in 
an optical system without any external feedback. Au- 
thors have predicted that in the scalar approximation, 
CP waves interacting in a NL Kerr medium character- 
ized by a noninstantaneous response, can undergo oscil- 
latory and chaotic temporal evolution, above a certain 
input intensity threshold. When the vector nature of 
light is included in the theoretical consideration it 
was found that the polarizations of CP light waves in an 
isotropic Kerr medium become temporally unstable, as 
the total intensity exceeds a certain threshold. Periodic 
and chaotic temporal behavior can occur in the output 
polarizations, as well as in the output intensities. Tem- 
poral instabilities in the polarizations of CP laser beams 
in atomic sodium vapor were investigated in p|. For in- 
tensities slightly above the instability threshold, the po- 
larizations fluctuate periodically. For higher intensities, 
the fluctuations become chaotic and the system evolves 
on a strange attractor. 

Continuous-wave and oscillatory transverse instabili- 
ties were predicted for CP waves in Kerr media for 
both the focusing and defocusing nonlinearities; neither 
cavity nor finite response time were required. Temporal 





FIG. 5: Counterpropagation of two beams in a 1 mm long 
crystal. The crystal is tilted by a = 10° with respect to the 
propagation direction, (a) shows the evolution of the forward 
beam (propagating from bottom to top); (b) shows the back- 
ward beam (propagating in the opposite direction), (c) shows 
the profile of the forward beam as it leaves the crystal. In (d) 
the dashed line shows the backward beam leaving the crystal 
after linear propagation, whereas the solid line shows it after 
nonlinear propagation. Adopted from [2^ ]. 



dynamics of the polarization state of CP waves in NL op- 
tical fibers was studied in [l^]. It was shown that in the 
presence of uniform twist, the dynamics may be reduced 
to an integrable chiral field representation. Dynamical 
instabilities of CP beams in a NL two-level system were 
investigated numerically in [10]. When the incident in- 
tensities are increased, this system becomes unstable and 
exhibits complex behavior, including quasi-periodic mo- 
tion and chaos. In [TtI ] the spatial polarization instability 
of two intense CP laser beams in an isotropic NL dielec- 
tric fiber was investigated experimentally. It was demon- 
strated that the distribution of polarization states along 
the fiber can be identified with a polarization domain 
wall soliton. 

Dynamical instabilities of CP self-trapped beams in 
PR media were reviewed in (28|. A route to chaos is 
described, including splitup instability, period doubling 
cascade, windows of intermittency, and fully developed 
chaos (Fig. [6]). An experimental method to stabilize 
unstable CP solitons using photonic lattices is developed 
by the same group; it is presented in the section dealing 
with the solitons in optical lattices. 



1 2 3 4 5 

FIG. 6: Bifurcation diagram, displaying transition to chaos 
in a ID model of CP self-trapped beams. Insets depict char- 
acteristic time dependence of the beams along the diagram, at 
one of the crystal faces. The steady beam (green) is the enter- 
ing beam, the unsteady beam (red) is the exiting beam. One 
can note steady-state CP soliton (upper left), single splitup 
transition (lower left), double splitup transition (upper mid- 
dle), period four oscillation (lower right), and chaotic response 
(upper right). Adopted from [2g| . 

III. TWO-DIMENSIONAL SYSTEMS 
A. Theoretical background 

Consideration of counterpropagation in two transverse 
dimensions in bulk media offers a more realistic and com- 
plete picture. Theoretical descriptions of CP self-trapped 
beams in 2D and time are provided in [25j and in [2?| . 
The derivation of equations is similar to the ID case [see 
Eqs. ([T])-(|3])]. the major difference being the appearance 
of the transverse Laplacian A x>a in 2D equations and 
different boundary conditions. However, the differences 
in physics and results are considerable, especially if one 
takes into account the anisotropic nature of the PR effect 
in 2D. We will confine our attention here to the isotropic 
approximation of the system in 2D. 

By assuming that the CP beams are incoherent (i.e. 
e = 0, making the E\ term disappear), the propagation 
equations are given by 

id z F + A x . y F = TE F , (6a) 
-id z B + A x<y B = TE B , (6b) 

and the temporal evolution of the system (or the time 
dependence of E ) is determined by Eq. (|3h.). 

The most important completely new effect in 2D, in the 
form of dynamical spontaneous symmetry breaking, was 
reported for the first time in [24j , where the counterprop- 
agation of two identical numerically calculated solitary 
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FIG. 7: CP fundamental beams for a medium of length 
L = 0.68-Ld displaying the splitup transition. The left col- 
umn shows the forward and the right column the backward 
beam (the arrows indicate the direction of propagation) . The 
top row is a snapshot after t — 25r. Both beams propagate 
through the medium as solitons. At t = 15GY (bottom row) 
the beams no longer propagate as steady solitons, but instead 
deviate on their way through the crystal. Reprinted from [2(| . 



beam profiles was investigated. For a value of T « 14, 
corresponding to a typical experimental situation, up to 
the length of the medium of 0.65Lrj no sign of instabil- 
ity is observed; the beams propagate as a CP soliton. 
However, at L = 0.68Lrj the solitary solution becomes 
unstable (Fig. [7]). At t = 25r the beams still propa- 
gate as solitons through their jointly induced waveguide, 
but the white noise included in the system excites an 
eigenmode that grows in time. At t — lOOr both beams 
start to deviate inside the medium from the straight ini- 
tial trajectories. Since the initial problem is rotationally 
symmetric, the direction into which the beams deviate 
is random in the isotropic approximation. The intensity 
distribution at t — 150t, presented in the bottom row 
of Fig. [71 shows a steady state of the system. This is 
an example of the splitup transition in 2D. The numeri- 
cal results show that the length of the medium and the 
power of the beams play an important role in the stabil- 
ity of CP solitons: decreasing/increasing the length (or 
the power) stabilizes/destabilizes the solitons. 

These results seem to contradict the results obtained 
for the solitons in copropagating geometry. Two mutu- 
ally incoherent solitons always attract each other, there- 
fore one would expect that the two CP beams always 
form a stable soliton. This is not the case. To find an 
explanation of the nature and the cause of the trans- 
verse splitup instability, the CP beams were considered 
as particles whose motion along the z axis is subject to 



the forces caused by the refractive index change in the 
medium. Because the medium is noninstantaneous, it 
was assumed that the motion of the " center of mass" of 
the beams is determined by the light distribution a time 
r ago. The second assumption was that the attractive 
force acting on the center of mass of each beam is propor- 
tional to the distance from the center of the waveguide 
induced in the medium by the beams. A simple har- 
monic oscillator-type theory of beam displacement that 
can account for the transverse shifts, derived in two in- 
dependent ways, was presented in [2J] and [27] ]. 



B. Pattern formation and linear stability analysis 

When excited beyond certain instability thresholds, 
very different physical systems display similar self- 
organized behavior that is described by the universal or- 
der parameter equations. A common necessary ingredi- 
ent is the MI of spatially uniform ground state, which 
leads to the spontaneous formation of extended periodic 
spatial structures. These patterns often exhibit simple 
geometric structure, such as rolls, rhombi, and hexagons. 
Linear stability analysis (LSA) provides a threshold for 
the static instability in such systems. 

NL optical materials are well suited for the observation 
of transverse Mis, especially in the CP geometry. The 
first complete CP pattern formation considerations and 
LSA in 2D was given in [2i| for the counterpropagation 
in a Kerr medium (in ID see [HI). ^ was demonstrated 
by a NL perturbation analysis that two very different 
pattern-forming modes coexist in this system. One is 
a hexagon-forming mode, and is dominant in the self- 
focusing media. The other is a roll-pattern mode, but 
it was found that rolls are unstable. Instead, square 
patterns emerge, and seem to be dominant in the self- 
defocusing media 1291 ■ When the two beams are slightly 
frequency-detuned [14j |, counterpropagation in PR two- 
wave mixing also gives rise to the transverse MI. The 
patterns that develop from the initial stage of MI are 
found to be predominantly rolls. For induced slight mis- 
alignment, full hexagonal patterns develop. 

The described phenomena are much more dependent 
on the geometry than on the particular form of the non- 
linearity. We will present here only one recent result 
concerning the transverse splitup instability of CP soli- 
tons (see Fig. [7]). Patterns developing in wider hyper- 
Gaussian CP beams will be covered in Sec. 6. In the 
standard MI theory one follows the dynamics of weak per- 
turbation to a wave and looks for instances of exponential 
growth of the perturbation. Such a growth promotes the 
amplification of sidebands and leads to the appearance 
of localized transverse structures. This approach is used 
much in the theory of transverse optical patterns [l5| . 
Here however, the whole object - a CP soliton - under- 
goes a sudden transverse shift to a new position. Using 
LSA, the splitup instability is explained as a first-order 
phase transition, caused by the spontaneous symmetry 
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breaking, and the threshold curve is determined [30(. 

One should note that LSA is more properly applied 
to very broad CP beams. In the case of splitup insta- 
bility the stability analysis is applied to a low-aspect- 
ratio geometry, and we are aware of its limited validity. 
It is known in many systems with dissipative feedback 
that the instability of solitons and pattern forming sys- 
tems follow different bifurcation routes. Here the insta- 
bility of propagating solitons and the pattern formation 
in wide CP beams (addressed in Sec. 6) are approxi- 
mately treated by the same LSA and the same threshold 
conditions. Qualitative agreement is found. 

One starts at the steady state plane-wave solution of 
the system of Eqs. ((3^) and (J6]): 




Fo(z) = F (Q)e- irE ° z , B (z) = B (L)e tFE ^- L \ (7) 

where E = -h/{\ + I ) and I = \F \ 2 + \B \ 2 . The 
primary threshold is determined by the linear instability 
of the steady state plane- wave field amplitudes Fo(z) and 
Bo(z), and the homogeneous part of the space charge field 
Eq. To perform LSA, a change of variables is made: 

F = F (l + f), B = B (l + b), E = E (l + e), (8) 

along with the change in the boundary conditions /(0) = 
b(L) = 0. Neglecting higher harmonics and terms 
quadratic in the perturbations /, b and e, and follow- 
ing the procedure described in Ref. [14], the threshold 
condition is obtained in a form: 



FIG. 8: Threshold curves obtained from Eq. (5). Reprinted 
from [Ij. 
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f^i sin *i sin* 2 „ ,^ 
l + cos* lC os* 2+ (^i + ^j L » =0 , (9) 

where *i = k 2 L, * 2 = \/fc 4 L 2 - 4ATk 2 L 2 , k being the 
transverse wavenumber. We chose |-Fo| 2 = I-6l| 2 j so that 
A = l^o| 2 /(l + 2|F | 2 ) 2 . This equation has the same 
form as the threshold condition in Ref. [2i|, except that 
the form and the meaning of variables and \&2 is dif- 
ferent. 

For each value of A there are two values of |F | 2 (or 
|i?i| 2 ) on the threshold curves represented in Fig. [8] 
For this reason we found it more convenient to plot the 
threshold intensity as a function of the square of the 
transverse wave vector (Fig. [9J; for each pair of values 
of r and L then one obtains different threshold curves. 
Also provided in Fig. [9] are the arrows which depict how 
much the CP solitons jump transversely in the k space 
in numerical simulations, after a splitup transition. The 
left end of an arrow points to the peak value of k 2 in the 
steady state, the right end points to the maximum value 
of the total transient change in k 2 . The end points are 
calculated by independent numerical runs of the full sim- 
ulations. For the given control parameters (r = 4 and 
L = 5Ljj) only single or double splitup transitions are 
observed. It is seen that the arrows provide a qualitative 



FIG. 9: Threshold intensity versus the square of the trans- 
verse wave vector k 2 , for F = 4 and L — 5Lo- Arrows cover 
the regions of jump in k 2 of the solitons in the inverse space, 
obtained numerically. Reprinted from [30t | - 

agreement with the form and the position of the lowest 
branch of the threshold curve, which signifies the first 
splitup transition. 



IV. EXPERIMENTAL DEMONSTRATIONS 

A. One-dimensional solitons 

The first experimental observation of spatial vector 
solitons in the counterpropagation geometry and for co- 
herent optical fields was reported by Cohen et al. [3l| . 
The experimental setup is shown in Fig. [11)1 An Ar + laser 
beam at 488 nm is split equally into two beams, 1 and 
2, that are focused to narrow stripes (15-mm FWHM) 
on the opposite faces, A and B, of an SBN:60 crystal 
(4.5 mm x 10 mm x 5 mm) in the configuration used ear- 
lier for the generation of PR screening solitons [32] . Two 
cameras image the two faces of the crystal. Importantly, 




FIG. 10: Experimental setup for the generation of CP soli- 
tons m. 



the light gathered by each camera consists of both the 
transmitted beam and the back-reflected beam. 

Figures [TITa.b) show the results of Ref. [3l[ for the 
image and intensity profiles taken by camera A at the in- 
put face B and at the output face A, respectively, when 
beam 1 is blocked and the nonlinearity is off. No soliton 
is formed. When the two beams propagate together with 
the nonlinearity (with an external voltage of 900 V), the 
beams mutually self-trap, as shown in (c). The FWHM 
of this combined beam is 15 /xm, equal to the FWHM of 
each of the input beams at both surfaces. Thus the com- 
bined wave, consisting of both CP beams, forms a vector 
soliton at the specific value of the nonlinearity, deter- 
mined by the applied field, the intensity ratio, and the 
crystal parameters (refractive index and the electro-optic 
coefficient). To exemplify the fact that the vector soli- 
ton is formed by both CP components, Cohen et al. [3l| 
blocked beam 1 and observed the output of beam 2 with- 
out changing the voltage. 

Further experimental studies of collisions between PR 
spatial solitons propagating in the opposite directions, 
performed by the same group [331 ] . demonstrated that 
each of the interacting solitons significantly affects the 
self-bending of the other, exhibiting effective attraction 
for one beam and repulsion for the other. In particular, 
Rotschild et al. [33| were able to switch between coherent 
and incoherent interactions by introducing a piezoelec- 
tric (PZ) mirror into the experimental setup, affecting 
one of the optical paths. Importantly, by varying the 
distance between the beams, the authors did not observe 
noticeable changes during the coherent collision between 
solitons. The coherent effects that occur during the colli- 
sion arise from the interference between the beams, trans- 
lated into a reflection grating. Rotschild et al. [33j tested 
the presence of such a reflection grating when the soli- 
tons are truly counterpropagating, by blocking one of the 
beams. The existence of this grating proved the occur- 
rence of a stable coherent interaction between the CP 
beams. When the PZ mirror is vibrating, the grating 
does not form; i.e., the soliton interaction is incoherent, 
and no reflection is observed from the grating. 




FIG. 11: Experimental images and intensity profiles taken 
by camera A (see Fig. I10p . Intensities of beam 2 in linear 
medium for (a) the input and (b) the output surfaces of the 
crystal when beam 1 is blocked (c) Total intensity of the CP 
vector solitons at the left face of the crystal, (d) Intensity at 
the left surface when beam 1 is blocked and the nonlinearity 
is on (adopted from Ref. [3ll]). 



B. Solitons in bulk media 

The first experimental study of CP solitons in a bulk 
medium was conducted by Jander et al. [34| . who also 
observed a dynamical instability in the interaction of CP 
self-trapped beams in a PR medium, predicted earlier 
in the theoretical modeling of the time-dependent beam 
dynamics (23|. Jander et al. [34} noticed that, while the 
interaction of copropagating spatial optical solitons ex- 
hibits only transient dynamics and eventually results in 
a final steady state, the CP geometry demonstrates a 
dynamical instability mediated by an intrinsic feedback. 
Experimental observations were found to be in qualita- 
tive agreement with the numerical simulations. 

Jander et al. [3_4j studied the dynamics of mutually in- 
coherent CP solitons in cerium-doped strontium barium 
niobate (Ce:SBN:60) crystal, using experimental setup 
shown in Fig. [T2j The crystal is biased by an external 
dc field along the transverse x direction, coinciding with 
the crystallographic c axis. Both beams are obtained 
from a single laser source and rendered mutually inco- 
herent by a mirror oscillating with a period significantly 
shorter than the relaxation time constant of the PR ma- 
terial. Propagating in different directions, both beams 
individually self-focus, and the nonlinearity is adjusted 
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FIG. 12: Experimental setup for the study of instabilities of 
CP solitons ;34j. Two beams are rendered mutually incoher- 
ent with an oscillating PZ mirror and focused on the opposite 
faces of a Ce:SBN60 crystal. Both crystal faces are imaged 
onto a CCD camera, allowing for synchronous observation of 
reflections of both the exit and input beams (Ms, mirrors; 
Ls, lenses; PH, pinhole; PBS, polarizing beam splitter; BS, 
beam splitter). Inset: CP soliton interaction in the numerical 
model (see details in Ref. |34|). 
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FIG. 13: Temporal plot of system dynamics, (a) Below 
threshold (Li=5 mm), the resulting stable and stationary 
state consists of two symmetrically overlapping solitons. (b) 
Above threshold (1/2=23 mm), irregular dynamics is ob- 
served 13411. 



such that each of the beams individually forms a spatial 
soliton. To demonstrate both above and below threshold 
behavior with a single crystal sample, Jander et al. (34j 
utilized two medium lengths by rotating the crystal about 
its c axis, thus yielding Li—5 mm and L2=23 mm. 

Initially, both beams are adjusted such that their in- 
puts and outputs overlap on both ends of the crystal, if 
propagating independently and in a steady state, includ- 
ing the shift through the beam bending. This configura- 



tion was chosen to minimize the possible effects of beam 
bending on the stability of a fully overlapping state in 
the form of a CP vector soliton, shown in Fig.[T3Ta). For 
comparison with the numerical simulations, experimen- 
tal data are reduced to one transverse dimension: The 
images obtained on the exit faces of the crystal are pro- 
jected onto the x axis. As these data are plotted over 
time, one gets a representation of the dynamics of the 
beam exiting the crystal face [as shown in Figs. [13ja,b)]. 
Changes parallel to the y axis are not represented, since 
most of the observable dynamics is confined to the x axis, 
owing to the significance of the c axis for the PR effect. 

For short propagation length (Li=5 mm), the output 
beams on both crystal surfaces initially shift their po- 
sitions, converging to an overlapping steady state, the 
vector soliton, see Fig. [T^a). In the case of a signifi- 
cantly longer medium (L2—23 mm) the beams initially 
self-focus separately [Fig. [T37 b)] and attract and over- 
lap. However, this state is unstable and yields to irregu- 
lar repetitions of repulsion and attraction. This process 
does not feature any visible periodicity and is observed 
for time spans that are orders of magnitude longer than 
the time constant of the system. 

More detailed experimental and numerical studies by 
Petrovic et al. [13] revealed that the CP incoherent beams 
can form 2D spatial solitons, but they also exhibit an in- 
teresting dynamical behavior in both dimensions. Stable 
solitons are readily observed over the 5 mm propagation 
distance, with an applied field of 1.3 kV/cm and the ini- 
tial beam peak intensity of about twice the background 
intensity. When the propagation distance is increased 
from 5 mm to 23 mm, for identical other conditions, 
a symmetry breaking transition from the stable over- 
lapping CP solitons to unstable transversely displaced 
solitons is observed. The beams still self-focus approx- 
imately into solitons, but they do not overlap anymore 
[Figs. rHT a.b)]. At the exit face most of the beam inten- 
sity is expelled to a transversely shifted position (about 
1 beam width), while a fraction of the beam remains 
guided by the other beam. This is another evidence of the 
splitup transition. At higher applied fields (stronger non- 
linearity) the beams start to move. The motion is such 
that the exiting beam rotates about or rapidly passes 
through the input beam, or dances irregularly around. 
No such long-lasting temporal changes are observed in 
the copropagation geometry. All these dynamical phases 
can qualitatively be reproduced by numerical simulation, 
as discussed above, in qualitative agreement with the ex- 
perimental results. 

The existence and stability of CP dipole-mode vec- 
tor solitons in a PR medium was studied experimentally 
by Schroder et al. 35], who also investigated the tran- 
sient dynamics. A dipole-mode vector soliton consists of 
two mutually incoherent beams: an optical dipole and a 
fundamental-mode beam. The individually propagating 
dipole does not form a spatial soliton, owing to repulsion 
of the dipole components. However, if a fundamental- 
mode beam that is incoherent to the dipole is launched 
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FIG. 14: Counterpropagating soliton after a splitup transi- 
tion: Forward propagating component in the steady state, (a) 
Exit face of the crystal, experimental, (b) The corresponding 
numerical simulation 13011 - 



between the dipole components, they will be trapped by 
means of incoherent attraction. 

Schroder et al. [35[ have proven the existence of a sta- 
ble CP dipole-mode vector soliton in a PR crystal. This 
vector soliton differs considerably in many aspects from 
its counterpart in the copropagating geometry. For ex- 
ample, the time scale of the transient dynamics is sig- 
nificantly longer. During the formation the beams split, 
and a trapped and an untrapped part of the dipole can 
be observed experimentally [351 ] . 



V. SOLITONS IN OPTICAL LATTICES 

Currently, we witness an explosion of interest in the 
localized CP beams in photonic crystals and optically 
induced photonic lattices pjj . Both ID and 2D geome- 
tries are extensively studied in a number of recent pub- 
lications, and many aspects of counterpropagation are 
addressed. An incomplete list includes Tamm oscilla- 
tions and localized surface modes at the interface between 
homogeneous medium and waveguide arrays; instabili- 
ties and stabilization influence of optical lattices; soliton 
propagation aided by reflection gratings; the behavior of 
CP vortices upon propagation in photonic lattices; Bloch 
oscillations and Zener tunneling; the questions of angular 
momentum transfer, conservation, and non-conservation 
in such systems; etc. In a few subsections, we will review 
some of these aspects. 



When the input power is increased, soliton instability 
sets in, and for a sufficiently high value the spontaneous 
symmetry breaking occurs, resulting in a discrete lateral 
displacement of the two solitons. A further increase in 
power leads to the temporally irregular behavior, with 
fast spatial fluctuations as compared to the buildup time. 
Numerical modeling [36j corroborated the existence of 
the three different regimes: stable propagation of vector 
solitons at low power, instability for intermediate power 
levels, leading to a transverse shift of the two discrete 
solitons, and an irregular dynamical behavior of the two 
beams at high input powers. 

The existence of Tamm oscillations at the interface be- 
tween a homogeneous medium and a ID NL waveguide 
array, with either cubic or saturable, self-focusing or self- 
defocusing nonlinearity, was demonstrated in [37[ . Light 
gets trapped in the vicinity of the edge of the array, ow- 
ing to the interplay between repulsion at the edge and 
Bragg reflection. Approximate analytical expressions for 
the repulsive potential are given for different types of 
nonlinearities. Tamm oscillations reduce to the surface 
Bloch oscillations, when the repulsive potential is a lin- 
early decreasing function of the distance from the edge 
of the waveguide array. 

The impact of an optically-induced photonic lattice 
on the dynamics of CP solitons in a biased PR crystal, 
as well as the stabilization of CP solitons by photonic 
lattices, were studied in [3^|. Numerical results there 
demonstrate that an optically-induced lattice of an ap- 
propriate period and strength can suppress or even com- 
pletely arrest the instability. It is found that CP soli- 
tons launched both on-site and off-site can be stabilized. 
The rate of decrease of temporal dynamics of CP solitons 
launched on-site strongly depends on the lattice strength 
and its period. In the case of small periodicity and high 
peak intensity of the lattice, spatiotemporal oscillations 
are observed, with characteristics dissimilar from those 
exhibited by CP solitons in bulk media. Experimental 
results [38[ demonstrate that, in most cases, soliton dy- 
namics, i.e. spatiotemporal oscillations, are suppressed 
with the increasing strength of the optical lattice. Also, 
the decrease in dynamics is noted experimentally for the 
ID lattice whose periodicity is comparable to the beam 
diameter. 



B. Solitons aided by reflection gratings 



A. One-dimensional lattices 

Interaction of CP discrete solitons in a nonlinear ID 
waveguide array was investigated experimentally and nu- 
merically in [361 ]. Solitons of equal input powers were 
launched in the same channel, but propagating in the 
opposite directions. Similar scenario to the CP soliton 
propagation in bulk is observed. For small input pow- 
ers the interaction between solitons is weak, and almost 
independent propagation in the same channel is seen. 



The existence of linear ID reflection grating in a 
NL optical medium implies periodicity along the optical 
propagation direction. In the presence of an intensity- 
dependent nonlinearity, the interference between the two 
CP beams produces a second longitudinal grating. Com- 
bining these effects with the standard diffractive broad- 
ening suppression, a twofold compensation mechanism 
can arise, allowing for the CP soliton formation. Soli- 
tons forming in a reflection grating in the presence 
of Kerr nonlinearity were investigated theoretically in 
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Bright and dark ID solitons counterpropagating along 
the reflection grating were analytically investigated in 
[39l | . It was shown that, depending on the Bragg match- 
ing between the light and the Fourier component of the 
grating, three different regimes of soliton existence arise. 
In the first regime, when the deviation from exact Bragg 
matching is small, both bright and dark solitons can ex- 
ist. The other two regimes occur for greater deviations 
from the exact Bragg matching. Deviations above Bragg 
matching allow only bright, and deviations below Bragg 
matching allow only dark solitons. In the two regimes the 
solitons are insensitive to the mutual phase difference. 

Transverse and soliton instabilities due to counterprop- 
agation through a reflection grating in Kerr media were 
considered in [Z(|. It was shown that the presence of 
the grating broadens and narrows the stability region of 
plane waves in focusing and defocusing media, respec- 
tively. Counterpropagation of spatial optical solitons in 
a linear reflection grating, in the presence of a longitu- 
dinally modulated Kerr nonlinearity, was investigated in 
[ill ] . The existence of symmetric soliton pairs supported 
by an effective Kerr-like nonlinearity is predicted ana- 
lytically. In addition, two families of solitons are intro- 
duced, in which the phase conjugation coupling exactly 
balances the Kerr holographic focusing. Properties of a 
general family of dark reflection solitons in defocusing 
Kerr media were considered in [12] • Two families of soli- 
tons counterpropagating along the grating direction in a 
modulated Kerr medium (asymmetrical one-solitons and 
coherent-incoherent two-solitons) were introduced in [43| . 



C. Two-dimensional lattices 

Self-trapped CP beams in fixed optical photonic lat- 
tices were for the first time investigated in j44[. When 
the propagation in photonic lattices is considered, the 
propagation equations are given by Eqs. © (assuming 
Eq — > E) , and Eq. ((3^) for the space charge field is mod- 
ified, to include the transverse intensity distribution of 
the optically induced lattice array I g : 



rd t E + E=- 



I + Ia 



1 + 1 + la 



(10) 



1. Launched vortices 

Rotational properties of mutually incoherent self- 
trapped vortex beams in optically induced fixed finite 
photonic lattices with a central defect were considered 
numerically in [45j. The defect is introduced either by 
omitting a lattice site in the center, or by defining a 
specific defect function. An interesting example of ro- 
tation of vortex filaments in the presence of a defect in 
the triangular photonic lattice is presented in Fig. 1151 
Although it looks as if the vortex filaments are rotating 
only within the defect in the center of the lattice, they 
also rotate away from the center, by tunneling between 
the lattice sites. We call this tunneling the nonlocal ro- 
tation, as opposed to the local rotation within the de- 
fect. The tunneling rotation is corroborated by the fact 
that, for all the cases presented in Fig. [T51 the angu- 
lar momentum of the vortex calculated over the whole 
lattice is considerably greater than the angular momen- 
tum of the vortex calculated only over the central part 
of the lattice. Nonlocal rotation in a periodic array, such 
as the triangular/hexagonal photonic lattice, can exists 
only for some values of the propagation distance L. Lat- 
tice supports nonlocal rotation only for some values of 
the propagation distance, with the "period" equal to 
(Fig. [15]). For the propagation distances between these 
values, chaotic behavior is observed, as well as the non- 
propagating modes. 

Rotating CP structures in the non-periodic but rota- 
tionally symmetric circular photonic lattices were dis- 
cussed also in [45[ . Results for the head-on counterpropa- 
gation of two centered vortices with the opposite topolog- 
ical charges in a circular photonic lattice, with a negative 
defect in the center, indicate rich dynamical behavior as 
a function of the control parameters T and L. For lower 
values of T or L stable structures are seen, in the form of 
well-preserved vortex core, centered at the defect, and fil- 
aments focused onto the neighboring lattice sites. Above 
this region stable rotating tripoles and quadrupoles ex- 
ist. For higher values of the parameters irregular rotat- 
ing structures and unstable structures (i.e. constantly 
changing structures of unrecognizable shape) are identi- 
fied. The rotating structures with filaments pinned to 
the lattice sites can exist only in the presence of lattice 
and have no analogs in the CP vortices propagating in 
bulk media. 



where / =|F| 2 + |i3| 2 is the total beam intensity. Spon- 
taneous symmetry breaking of the head-on propagating 
Gaussian beams is observed, as well as discrete diffrac- 
tion and the formation of discrete CP vector solitons. In 
the case of launched vortices, beam filamentation is re- 
ported, and subsequent pinning of filaments to the lattice 
sites J44|. Dynamical properties of mutually incoherent 
self-trapped CP beams in optically induced photonic lat- 
tices, for different incident beam structures and different 
lattice configurations, were investigated numerically in a 
number of papers 0, HI S3, H \M S3, EH ■ 



2. Gaussian beams 

Time-dependent rotation of CP mutually incoherent 
self-trapped Gaussian beams in periodic optically in- 
duced fixed photonic lattices was investigated numeri- 
cally in [H, |43] . In these papers lattice arrays with the 
square or triangular arrangements of beams were con- 
sidered, with the central lattice beam absent. Head-on 
CP Gaussian beams were launched into the center of the 
lattice, parallel to the lattice beams. For both photonic 




FIG. 15: Local and nonlocal rotation of vortices in triangu- 
lar lattice for various propagation distances. For propaga- 
tion distances between these values chaotic behavior is ob- 
served. Input vortices have the opposite topological charge 
±1. Intensity distribution of the forward field at its out- 
put face is presented for parameters: F = 17, lattice spacing 
d = 28fim, FWHM of lattice beams 12.7^m, input FWHM 
of vortices 26.2fim, maximum lattice intensity I g — 5Id, 
\F \ 2 = \B L \ 2 = 5I d . Reprinted from g|. 

lattices the periodic rotation was found (see Fig. [TH] for 
triangular lattice) in a very narrow region of control pa- 
rameters. Each Gaussian beam collapses to a displaced 
soliton-like beam, and after transient dynamics starts to 
rotate indefinitely. Since for the parameters of such sta- 
ble periodic solutions there exist no stable steady states, 
and since in numerics Eq. (|10p becomes equivalent to 
the scalar nonlinear delay-differential equation, this phe- 
nomenon is recognized as a supercritical Hopf bifurca- 
tion. The central parts of Gaussians rotate regularly in 
the center of the lattice, owing to the defect, and along 
the whole crystal. Filaments away from the center rotate 
with the same frequency about the symmetry axis of the 
lattice, by tunneling between the lattice sites, but only 
close to the exit face of the crystal. 

Gaussian-induced rotating structures present soliton 
solutions, because they preserve shape along the main 
symmetry direction during the rotation. The physical 
origin behind the nonlocal rotation is incoherent inter- 
action and spontaneous symmetry breaking, while the 
rotation is realized through the tunneling. Observed ro- 
tating structures are stable in the presence of up to 5% 
noise added to the input beam intensity and phase. Spon- 
taneous symmetry breaking via noise determines the di- 
rection of rotation, both directions occurring with 50% 
probability. For the same control parameters, CP Gaus- 
sian beams show very irregular dynamical behavior in the 
absence of lattice, but very stable propagation is found 
in the presence of lattices without defects. 



3. Stable rotating solitons 

For geometries and parameters which allow stable ro- 
tation, the existence of solitonic 2D solutions was investi- 
gated by considering Eqs. (O and (JTQJ) in the steady state 
|46| . Due to their symmetry, the above equations suggest 
the existence of a fundamental 2D soliton solution of the 
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FIG. 16: Gaussian-induced rotation for the triangular pho- 
tonic lattice: intensity distribution of backward beam at its 
exit face of the crystal, presented at different times. The 
second row shows isosurface plots of a rotating structure at 
characteristic times. Parameters: F = 25, input FWHM of 
CP Gaussian beams ll^m, beam power 3.80, propagation 
distance L = 2Ln — 8mm, input beam intensity |-Fb| 2 = 
\B L \ 2 = lid, lattice spacing d = 28/im, FWHM of lattice 
beams 12.7^tm, and maximum lattice intensity 7o = 20/d . 
Adopted from [461 ]. 

form: 

F = u(x, y) cos(6>)e' iAlz , B = u(x, y) sin(6>)e~ vz , (11) 

where fj, is the propagation constant and 9 is an arbi- 
trary projection angle. An appropriate choice for the CP 
beams is 9 — 7r/4; the same analysis can be applied to 
the copropagating geometry, but with the choice 9 = 0. 
When this solution is substituted in Eqs. ©, they both 
transform into one, degenerate equation: 

l?;l 2 4- I 

(M + Au + Tu I 7 g =0. (12) 

1 + M 2 + Ig 

The soliton solutions of Eq. (|T2"j) can be found using the 
modified Petviashvili iteration method. Because of the 
CP geometry, these solitonic solutions are stable only up 
to some critical value of the propagation distance. Gaus- 
sian input beams and the same parameters as in the full 
numerical simulations are used in search of the stable 
soliton solutions. The propagation constant [i is varied, 
in order to find the beam power (P = J J \u\ 2 dxdy) cor- 
responding to the stable rotating structures. Figure [17] 
depicts the power diagram, together with the character- 
istic soliton solutions for the case of triangular photonic 
lattice; the filled circles represent the characteristic types 
of symmetric discrete soliton solutions. By increasing the 
propagation constant [i these solutions become more lo- 
calized and asymmetric. Only for the beam powers cor- 
responding to the less localized symmetric soliton solu- 
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FIG. 17: Power diagram of soliton solutions for the triangular 
photonic lattice. Various symbols indicate different kinds of 
soliton solutions, characterized by the corresponding profiles. 
Reprinted from [461 ]. 



tions, and for the lower values of /z, can one find Gaussian 
induced rotation in numerical simulations. 



4- Angular momentum transfer 

The transfer of orbital angular momentum (AM) from 
vortex beams to an optically induced photonic lattice was 
demonstrated numerically in [5l|. An optically induced 
photonic lattice is, in fact, a complex propagating laser 
beam, which means that the propagation equations for 
the total system of interacting incoherent CP beams in 
the computational domain should be of the form: 



id z F = -AF + TEF, id z G f = -AG/ + TEG f , (13a) 
id z B = AB — TEB, id z G b = AG b - r£G fc ,(13b) 

where G/ and Gb are the envelopes of the forward and 
backward propagating lattice beams, and / = \F\ 2 + \B\ 2 
and I g = |G/| 5 



|Gb| are the corresponding beam in- 
tensities; the temporal evolution of the space charge field 
is given by Eq. ([IT))) , as before. 

Numerical results for the transfer of AM in the inter- 
acting CP beams are presented in Fig. [TH1 It was found 
that the transfer of orbital AM is minimal in the inter- 
acting CP lattices, and that the total AM - meaning the 
sum Lz(F) + Lz(Gf) + Lz(B) + Lz(Gb) of all momenta 
along the propagation z axis - is not conserved in this 
case. More precisely, the difference of AM of CP beams 
Lz(F) + Lz(G f) — Lz(B) — Lz(Gb) is conserved, whereas 
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FIG. 18: Transfer of angular momentum in the interacting 
CP beams, (a) Forward lattice at the input and output faces, 
and the vortex at the output face, (b) Different normalized 
angular momenta. Parameters: r = 3, L = 2.5Ld = 10 
mm, lattice spacing 28/im, FWHM of input vortices 24.6/im, 
FWHM of input lattice beams 9/im, maximum input intensity 
\Fo\ = \Bl\ = 5Id, and the maximum input lattice intensity 
\G f \ 2 (z = 0) = \G b \ 2 (z = L) = 20I d . Adopted from [EH]. 



their sum is not. Even though any real optically induced 
photonic lattice is an interacting beam, in a number of 
papers it is treated approximately as a fixed lattice. It 
was found in [5l| that the transfer of orbital AM can be 
substantial in the fixed periodic lattices, and that there 
is always a considerable loss of AM. Only in the fixed ra- 
dially periodic lattices there is no problem with the con- 
servation of AM of propagating light - it is a conserved 
quantity there. 

Different behavior noted in the interacting and fixed 
lattices can rigorously be explained [5l|. By using the 
standard definition for the z component of the orbital 
AM, Lz(F) = -i/J dxdyF*(x,y)(xd y - yd x )F(x,y) + 
c.c, the derivative of the difference of AM for the as- 
sumed CP geometry of propagation in the steady state 
is given by: 

dLz(F) dLz(G f ) dLz(B) dLz(G b ) 



dz 



dz 



dz 



dz 
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r 2 * d(\F\ 2 



|G b | 2 ) 



dip 



pd P T[ln(l + I + I g )-I-I g ] 



271 

) 





(14) 



where p and ip are the cylindrical coordinates. The differ- 
ence of AM is conserved in the case of interacting CP lat- 
tices, because the space charge field E is then an explicit 
function of / + I g , the integration in ip is over a perfect 
derivative, and the integral in Eq. (jT4j) is 0. However, 
for the fixed periodic lattices the terms involving G/ and 
Gf, are absent, while E still contains I + I g , and the in- 



tegration in Eq. (| 14[) is not over a perfect derivative, so 
the integral does not vanish. Therefore, the difference 
in AM is then not conserved. (Note that for the fixed 
radially periodic circular lattices the difference in AM is 
conserved.) The same Eq. (fl"4]) proves the nonconser- 
vation of the sum of AM in the general CP case: the 
integral then contains the difference of intensities, while 
E still contains the sum, so it can not be equal to 0. 



D. Experimental demonstrations 

A simple realization of a periodic NL medium is the 
one-dimensional system, in which the NL arrays con- 
sist of parallel, weakly coupled waveguides. In these ar- 
rays, discrete soliton interaction has been investigated 
for parallel beams, showing soliton attraction, repulsion, 
oscillatory behavior of the two beams, and soliton fu- 
sion [5^, HH . Similar to what has been observed for spa- 
tial CP solitons in bulk media, both instability of the 
interacting discrete solitons, leading to discrete spatial 
shifts, and irregular dynamical behavior for high non- 
linearities can be observed in photonic lattices. Experi- 
mental results of Smirnov et al. [3a | reveal the existence 
of three regimes, namely, the stable propagation of vec- 
tor solitons, an instability regime leading to discrete dis- 
placements of solitons, and an irregular dynamics regime. 

Figure [TU] presents numerical calculations which sum- 
marize these three regimes. For a small intensity ratio the 
two solitons propagate stably with only weak interaction, 
whereas for a higher ratio instability grows and the soli- 
ton formation is partly suppressed. For even higher inten- 
sity ratio the two formed discrete solitons are displaced 
by one channel to the left. If the intensity ratio is fur- 
ther increased, no steady-state solution can be obtained 
anymore: The output intensity on both faces starts to 
fluctuate rapidly over recording times, similar to the re- 
sults described in Ref. [25j. In experiment, the interac- 
tion of discrete CP solitons has been studied by Smirnov 
et al. [361 ] in LiNbC>3 waveguide array. For small input 
powers or intensity ratios, respectively, the interaction is 
weak and almost independent propagation of the two dis- 
crete solitons in the same channel is achieved. When the 
input power (intensity ratio) is increased, soliton insta- 
bility occurs, and for sufficiently high values spontaneous 
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FIG. 19: Numerical simulation of the interaction of discrete 
CP solitons (left-hand side, beam 1; right-hand side, beam 2) 
for three different intensity ratios [36f ] . 



symmetry breaking with discrete lateral displacement of 
the two solitons is observed. 

These observations led to the extension of efforts on 
stabilization of CP solitons in 2D volume systems b y th e 
use of ID and 2D photonic lattices (Koke et al. [3a|). 
They investigated the dependence of the instability dy- 
namics on the period and amplitude of the lattice, and 
presented experimental verification for the dynamic sta- 
bilization of the bi-directional soliton states. Fig. [20] 
displays the arrest of instabilities of 2D CP solitons in 
ID photonic lattice. It is evident that the ST oscillations 
are suppressed by increasing the strength of the optical 
lattice. This evolution is accompanied by an increased 
trapping of light in the neighboring lattice sites, as ob- 
served in the case of ID CP solitons. In contrast to the 
ID case, an asymmetry of trapping in the lattice chan- 
nels is due to the beam-bending effect. Owing to the 
change in the refractive index, the self-bent soliton gets 
partially reflected when it passes a lattice site, and the 
reflected light travels along the waveguide written by the 
lattice wave. Moreover, the weaker oscillations in the 
y-direction are suppressed as the lattice peak intensity 
is increased. It is worth noting that the experimental 
stabilization of CP solitons has been achieved with lat- 
tice strengths much lower than that found in numerical 
simulations. 

In experiments, Koke et al. [HI used a 2D square lat- 
tice, optically-induced in a PR crystal. They varied the 
lattice period and power, and monitored the positions of 
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FIG. 20: Temporal evolution of the intensity distribution of 
2D CP solitons in optical lattices, projected onto the ir-axis 
(parallel to the c-axis) for different lattice powers of (a) 250 
fiW, (b) 1.0 mW, (c) 1.5 mW, (d) 2.0 mW, and (e) 2.5 mW. 
The faint horizontal lines at x = mark the reflected beam 
at this crystal face, which acts as a reference. The results 
for the other crystal face show a similar behavior. Reprinted 
from HI. 
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FIG. 21: Stabilization of the instability in two transverse di- 
mensions: (a) two-dimensional square lattice of 6 /jm period, 
(b) Input intensity distribution for the forward (F) and the 
backward (B) propagating beams, respectively. (c,d) Digi- 
tally combined beam profiles at the left and right faces of the 
crystal, when each soliton propagates independently. (e,f) 
Stabilization by the lattice on both faces of the crystal, re- 
spectively [38l ]. 



beam self-focusing and self-bending. At bias electric field 
of 2 kV/cm and at powers of 1 mW each beam forms a 
spatial soliton, where the soliton size is equal to the input 
beam size. In Fig.[2TTc.d) we show the digitally combined 
input and output of each beam, as they would propagate 
without interaction inside the crystal. When both beams 
co-propagate, they start to interact. After the initial 
attraction, the beams exhibit oscillatory dynamics. 



VI. GENERALIZATIONS 



the beams at both faces of the crystal. The inputs of the 
forward and backward propagating beams in their exper- 
iments have the size of 19 /im and 18 ^tm, respectively 
[Fig. [2TTb')]. In 10 mm long crystal this corresponds to 
approximately 5 diffraction lengths of linear propagation. 
For a small lattice period (3 ^m) the potential induced 
by the lattice was too weak to arrest the instability of 
the CP beams. With the increased lattice period (of 6, 
9, and 12 /im) the instability was practically removed for 
a certain range of lattice strength. The large lattice peri- 
ods, however, strongly reduce the mobility of the beams, 
as each beam can be fully trapped at a single lattice site. 
Such trapping imposes a constraint on the formation of 
bi-directional waveguides, which becomes sensitive on the 
initial alignment of the beams. Thus, beams propagating 
in different directions inside the crystal will not attract, 
as their intensity overlap will be reduced by the trapping 
on different lattice sites. 

Without the lattice, both beams overlap weakly and 
their individual propagation is strongly affected by the 



A. Localized multipole beams 

Dynamical behavior of mutually incoherent CP mul- 
tipole vector solitons in an SBN:60Ce PR crystal was 
investigated in [54| , The dipole-dipole interaction from 
that paper is reproduced in Fig. [2H In the case of dipole- 
dipole CP solitons, two identical dipoles with the compo- 
nents out of phase were countcrpropagated head-on. The 
dipoles were aligned perpendicular to the external field, 
which points in the horizontal direction. A transverse 
split-up occurred, the direction of the split-up is prefer- 
entially along the direction of the external field, and it 
also depends on the added noise [Fig. [22T b)] ■ Only in the 
case when some noise is added to one of the beams was it 
possible to observe skewed split-ups, in better agreement 
with the experiment. For the case with no noise [Fig. 
[227 c)]. in the beginning oscillations were noticed along 
the y axis, and after a short time these oscillations were 
damped. Compared to the single CP soliton cases, the 
cases involving dipoles are more stable and the transient 
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FIG. 22: Dipole-dipole interaction, (a) Experiment: the for- 
ward beam (upper) and backward beam (lower), at the exit 
face of the crystal. The corresponding numerical simulations 
of the backward beam: (b) with an extra noise of 5% added 
to the input beam intensity, (c) without noise. Parameters: 
maximum input intensity |-Fb| 2 = \Bl\ 2 = 1.3/d, T = 7.17, 
L = 5.75L_d = 23 mm, initial beam widths (FWHM) 20/xm, 
and the initial distance between the dipole partners 40/im. 
Reprinted from [541 ]. 



dynamics last shorter. 

The development of higher-order multipole structures 
and patterns in CP beams in saturable Kerr-like me- 
dia was investigated in [5oT ]. A systematic numerical 
study was carried out, by varying the width of beams. 
The results of LSA, concerning the instability of plane 
waves, were compared with the numerical results con- 
cerning broad hyper-Gaussian beams (used as inputs in 
simulations) whose width was varied. Qualitative agree- 
ment was found, due to the similarity between the plane 
wave and the flat-top hyper-Gaussian beam profile. We 
should again stress the fact that the splitup transitions 
do not appear to be of this common type of MI. The 
solitons themselves could be considered as related to the 
filaments of MI, and, as such, should be stable against 
the same kind of MI. Nonetheless, it is still of interest 
to explore the cross-over region by increasing the size 
of the solitonic beams, until they display MI. A smooth 
transition from the soliton splitup instabilities of narrow 
beams to the pattern-forming transition of broad beams 
is observed. 

An interesting consequence of the finite size effects is 
the appearance of the circular saw instability, presented 
in Fig. 1231 It appears in the form of circular saw-like ro- 
tating blades, visible in an intermediate region of beam 
widths, and it is caused by the MI at the edge of the 
beam profile. It happens very close to the absolute min- 
imum of the control parameter ATL = 7r/4 and is very 
robust. The rotation of optical beams along the propaga- 
tion direction in saturable Kerr-like media generally ap- 
pears through a bifurcation in the spatial domain. Here 
the bifurcation from hyper-Gaussian beams into rotating 
structures is due to a spatial symmetry breaking associ- 
ated with a Hopf bifurcation in the time domain. 

For a higher value of the control parameter ATL = 2, 
where one cannot expect that LSA is applicable (see Fig. 
[8]) , a more complex behavior in the form of higher-order 
multipole structures was found (Fig. |2~4")) . For the nar- 
row width of incident beams, FWHM = 20/im, a rotat- 
ing displaced soliton was seen at the exiting faces of the 
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FIG. 23: Circular saw instability of the backward beam, pre- 
sented in the direct and in the inverse space (the second col- 
umn), for ^4ri/ = 0.8 (close to the absolute minimum of the 
threshold curve). Transverse intensity distributions of the 
backward beam at the exit face are presented at different 
times, for three values of FWHM, recorded at the left edge 
of each of the rows. Circular saw- like rotating blades be- 
come visible in a region of FWHM, after a long transient 
development. Parameters: |F(0)| 2 = \B(L)\ 2 = 2, T = 6.68, 
L — 1.5-Ld. The size of the transverse window in the direct 
space is 400 x 400/im. Adopted from [55|. 
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FIG. 24: Different dynamical multipole structures as the 
FWHM of the backward beam is increased, presented at two 
instances, for the control parameter ATL = 2. The fig- 
ure should be viewed as a two-column picture. Parameters: 
|F(0)| 2 = \B(L)\ 2 = 3, T = 10.9, L = 3L D . The size of the 
transverse window is 360 x 360/im. Adopted from [551 ] - 



crystal, after a split-up transition. At a larger width 
(FWHM = 100/im) a hexagonal structure was observed 
in the beginning, which was followed by a regular rotation 
of filaments. Since for the parameters of such stable peri- 
odic solution there exist no stable steady state, and since 
numerically Eq. ([3H) is equivalent to a scalar nonlin- 
ear delay-differential equation, this phenomenon is recog- 
nized as the supercritical Hopf bifurcation. For FWHM 
= 140^m, steady octagonal structure appeared. For the 
next FWHM, after a set of various regular patterns, ir- 
regular structures took place. It should be mentioned 
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FIG. 25: Iso-surface plots of two CP vortices, with charges 
±1. Upon collision the vortices break into three beamlets, 
which co-rotate continually in the sense indicated by the ar- 
rows. Adopted from [241 ]. 



that all of these structures appear after a prolonged tem- 
poral development, when the secondary instabilities set 
in and the interaction of NL modes takes place. 



B. Vortex beams 

Colliding vortices with opposite topological charges 
were considered in [25|. In this case the beams break 
up, generally into more than two filaments. In Fig. 
ESI the case of two CP counter-rotating vortices is dis- 
played, which break into three beamlets with phase shifts 
of 2"7r/3. After a while the beams form stable rotating 
structures that do not change in time. When viewed 
in their exit faces, the beams form true rotating pro- 
pellers. Such stable rotating state is interesting for devel- 
oping transverse Mis over a fraction of diffraction length, 
even though it is generated in an isotropic model. Previ- 
ously observed isotropic vortex vector solitons, with co- 
propagating components, tended to propagate for tens of 
diffraction lengths before developing Mis. 

Optical CP vortices in PR crystals were investigated 
numerically in [56]. A general conclusion of numerical 
studies there was that the CP vortices in a PR medium 
cannot form stable CP vortex (i.e. ring-like) structures, 
propagating indefinitely. For smaller values of T or the 
propagation distance L stable CP vortices were observed. 
Nevertheless, when they break, they form very differ- 
ent stable filamented structures in propagating over fi- 
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FIG. 26: Typical behavior of CP vortices in the parameter 
plane. The input vortices have the same topological charge 
+1, and maximum input intensities |-Fb| 2 = \Bl\ 2 = 3. Insets 
list the possible outcomes from vortex collisions. Reprinted 
from [561 ] . 



nite distances, corresponding to typical PR crystal thick- 
nesses, which are of the order of few Lr>. Numerical 
studies showed that the CP vortices with the same topo- 
logical charge tend to form standing waves, whereas the 
vortices with the opposite charges tend to form rotating 
structures. Some typical examples of collisions between 
single head-on input vortices with the same topological 
charge +1 are shown in Fig. fSHl which represents the 
phase diagram in the plane of control parameters. One 
can notice in the figure a narrow threshold region which 
separates the stable vortices from the other structures. 
The shape of the threshold region follows the general 
TL = const, form. Above this region stable dipoles, 
tripoles and quadrupoles are seen, in the form of standing 
waves. For higher values of the parameters, the following 
quasi-stable situations are identified: the transformation 
of a quasi-stable quadrupole into a stable tripole, several 
transformations of quadrupoles into quadrupoles, and a 
stable rotating dipole. Above the quasi-stable region, CP 
vortices produced unstable structures. 



C. Counterpropagating beams in liquid crystals 

Nematic liquid crystals (NLC) exhibit huge optical 
nonlinearities, owing to large refractive index anisotropy, 
coupled with the optically-induced collective molecular 
reorientation. They behave in a fluid-like fashion, but 
display a long-range order that is characteristic of crys- 
tals. Thanks to the optically nonlinear, saturable, non- 
local and nonresonant response, NLC have been the sub- 
ject of considerable study in recent years. The behavior 
of CP self- focused beams in bulk NLC, both in time and 
in three spatial dimensions was investigated in [5^ . [5c| , 
using an appropriately developed theoretical model and 
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a numerical procedure based on the beam propagation 
method. 

The evolution of slowly- varying beam envelopes F and 
B, linearly polarized along x axis and propagating along 
z axis in a NLC cell, is described by the following paraxial 
wave equations: 



3F 

2ik — + AF + /c 2 e a [sin 2 6 - sin 2 (<9 rcst )]F = , (15a) 
oz 

3B 

-2ik — + AB + fc 2 e a [sin 2 9 - sin 2 ((9 rcst )]B = , (15b) 

where F and B are the forward and backward propagat- 
ing beam envelopes, k = k^no is the wave vector in the 
medium, e a = n 2 — n 2 , is the birefringence of the medium, 
and A is transverse Laplacian. The rest distribution an- 
gle 9 Tes t m the presence of a low-frequency electric field 
is modeled by: 



0re S t{z,v) = e o (v) + [e ia -e o (v)] 



[exp(— z/z) + exp(— (L — z)/z)}, 



(16) 



with 9q(V)) being the orientation distribution due to the 
applied voltage far from the input interface. 9 m is the 
director orientation at the boundaries z — and z = L, 
where L is the propagation distance and z is the relax- 
ation distance. The temporal evolution of the angle of 
reorientation is given by the diffusion equation: 



09 

7-^- = KA X v t 

I gt x,y 



-e e a sin(20)[|F| 2 



\B\ 



(17) 



where 7 is the viscous coefficient and K is Frank's elastic 
constant. Here 9 is the overall tilt angle, owing to both 
the light and the voltage influence. 

It was found numerically [53, [58[ that the stable vec- 
tor solitons can only exist in a narrow threshold region 
of control parameters. Bellow this region the beams 
diffract, above they self- focus into a series of focal spots. 
Spatiotemporal instabilities were observed as the input 
intensity, the propagation distance, and the birefringence 
were increased. The effect of the input intensity variation 
on the CP Gaussian beam propagation is presented in 
Fig. For smaller intensities [Fig. 127(a)] self-focusing 
is too weak to keep the beam tightly focused, so it can 
not pass through unchanged, as a spatial soliton. By 
increasing the beam intensity [Fig. 127( b)] at one point 
stable soliton propagation is achieved. For still higher 
intensities transverse motion of the beam is observed, in 
the form of one [Fig. [2"7T c)], or two consecutive jumps, 
resembling beam undulations. For further increase of 
the intensity unstable dynamical behavior of beams [Fig. 
l27Td)] is seen. 

An interesting experimental account on the interaction 
of CP solitons in a NLC E7 cell is provided in Ex- 
periments are performed to estimate the nonlocality of 
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FIG. 27: Beam propagation, shown for one beam in the (y, z) 



plane, for different input intensities: (a) I 
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(b) I = 7 x 10 9 V 2 /m 2 , (c) I = 8 x 10 9 V 2 /m 2 , and (d) I = 
9 x 10 10 V 2 /m 2 . In (c) the beam intensity is also shown in the 
(x,y) output plane. Parameters: input beam width (FWHM) 
4/im, L — 0.5 mm, and e a = 0.5. Adopted from [BH ]. 
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FIG. 28: Curved beam resulting from the fusion of two nonlo- 
cal CP solitons. The fibres are artificially marked, the dashed 
lines indicate the initial direction of the solitons. Reprinted 
from HI. 
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the reorientational nonlinearity in thick samples. The 
attraction of spatial optical solitons counterpropagating 
parallel to each other and at different small distances 
is displayed. An experimental method to estimate the 
width of the refractive index profile in a NLC sample ex- 
cited by a narrow laser beam is developed. It is shown 
that the width of the index profile can nicely be fitted 
by a Lorentzian curve. A rare experimental picture of 
a stable CP soliton pair, launched with a considerable 
transverse displacement in a NLC, is presented in Fig. 

m 

VII. CONCLUSIONS 

We have summarized recent developments in the 
physics of CP optical beams and spatial solitons, propa- 
gating in NL media. We have analyzed the formation of 
various stationary modes, as well as spatiotemporal insta- 
bilities of CP beams. We have employed several models 
for describing the evolution and interactions of optical 



beams and spatial solitons that propagate in opposite di- 
rections, but the majority of the results are presented for 
the model of saturable PR nonlinearity. We have dis- 
cussed the recent experimental observations of the coun- 
tcrpropagation effects and instabilities in waveguides and 
bulk geometries, as well as for one- and two-dimensional 
photonic lattices. We have also discussed several gen- 
eralizations of this concept, including the CP beams of 
complex structures, such as multipole beams and optical 
vortices, as well as counterpropagation in other media, 
such as photonic and nematic liquid crystals. 
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